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Abstract 

We consider nearly-flat-band Hubbard models of a ferromagnet, that is the models 
that are weak perturbations of those flat-band Hubbard models whose ground state is 
ferromagnetic for any nonzero strength U of the Hubbard repulsion. In contrast to the 
flat-band case, in the nearly-flat-band case the ground state, being paramagnetic for 
[/ in a vicinity of zero, turns into a ferromagnetic one only if U exceeds some nonzero 
threshold value Uth- We address the question whether Uth of the considered models 
is in a physical range, therefore we attempt at obtaining possibly good estimates of 
the threshold value Uth- A rigorous method proposed by Tasaki is extended and the 
resulting estimates are compared with small-system, finite-size scaling results obtained 
for open- and periodic-boundary conditions. Contrary to suggestions in literature, we 
find the latter conditions particularly useful for our task. 



1 Introduction 

The origin of ferromagnetism in itinerant electron systems, that is of the existence of a net 
macroscopic magnetic moment, without any external magnetic field, is an old problem set 
up already by founding fathers of quantum mechanics. Following almost equally old hints by 
Heisenberg [Ij, it is the Fermi statistics and the Coulomb interaction between the electrons 
that are responsible for this phenomenon. In their quest for ferromagnetism in itinerant 
electron systems, Mielke and Tasaki (see [2j and references quoted there) have analyzed 
the paradigmatic Hubbard models that describe strongly-correlated electron systems. The 
general Hamiltonian of those models reads: 

H = ClC,,„ + ^uJ2n^,ani,.., (1) 

where the sums are over all the sites i,j of the underlying lattice, and over projections of the 
electron spin a on some axis. The operators Cj^, Ci^„, stand for the creation and annihilation 
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operator, respectively, of an electron whose spin projection is a, being in a state belonging 
to an orthonormal basis of the single-particle state space, that is localized at (and therefore 
labeled by) lattice sites i. The coefficients tjj - the matrix elements of a single-particle 
Hamiltonian between states \i) and \j) - give the hopping intensities (if i ^ j) and on-site 
external potentials (if i = j). The term proportional to f/ > is a two-body interaction that 
represents a strongly screened Coulomb repulsion; rii^a- = CjaC'i,a- By varying the underlying 
lattices, the matrices tij, and the number of electrons, A^e, we obtain different systems of 
itinerant electrons in tight-binding approximation. From the physical point of view, the 
Hamiltonians ([T]) describe extremely simplified systems of itinerant electrons. However, 
these models are not too simple to be relevant for the question of ferromagnetism; the 
crucial property, that is the invariance with respect to rotations in the total-spin space 
{SU{2) invariance in the total-spin space), holds true. 

Apparently, the two "ingredients" of ferromagnetism, suggested by Heisenberg, are built 
in these models. For suitable choices of the underlying lattice, the matrix elements tj j, and 
the electron number N^., Mielke and Tasaki have succeeded in proving [2j, for the first time, 
that the ground state is ferromagnetic, that is the total spin, 5*^*, is proportional to the 
electron number Ne, for any U > . Moreover, that ground state is a saturated ferromagnet 
(that is, for given Ng, Stot is maximal). 

Their results were the first ones providing us with some fine insight into the nature of 
itinerant ferromagnetism, with some sufficient conditions for ferromagnetism. Those results 
go far beyond the mean-field theory that leads to the qualitative Stoner criterium a kind 
of necessary condition for ferromagnetism. In their analysis, besides the Fermi statistics and 
Coulomb interactions, another factor played a crucial role - the topology of the underlying 
lattice. Without going into details, this topology should enable one to adjust the hopping 
intensities and external potentials in such a way that the lowest band in the single-particle 
spectrum is non-dispersive, or briefiy flat, and there exists a basis of the fiat-band eigen- 
subspace, which consists of localized (i.e with bounded support) eigenstates with suitable 
overlap and connectivity properties. 

Now, it has to be mentioned that the theoretical limit of fiat band results in a nonphysical 
system that cannot be realized in any experiment. This is because such an electron system 
violates a number of physical principles. For instance, in a range of number densities, the 
residual entropy is finite (see [3] for examples where it was calculated exactly), and there is 
no Fermi surface. Moreover, the density of states of fiat-band systems is singular and there is 
no competition between the kinetic energy and the Coulomb interaction energy. Therefore, 
as far as the physical phenomenon, the itinerant ferromagnetism, is concerned, the described 
models would amount for not more than toy models, if not the subsequent work of Tasaki 
([5] and references quoted there). 

Tasaki has been able to prove that the flat-band ferromagnetism of specific Hubbard 
models is robust against perturbations that broaden the fiat band. At least at some cases, 
the perturbed fiat-band models, the so called nearly-flat-band models, have also ferromagnetic 
ground states but only if the value of the Hubbard repulsion U is greater than some threshold 
value. For sufficiently large systems and narrow broadened fiat bands that threshold value 
does not depend on the kind of perturbation and boundary conditions applied to the system, 
and is, therefore, an intrinsic characteristic of nearly-fiat-band electron systems; in the sequel 
denoted Uth- This quantity appears to be the greater the stronger the perturbation is, for 
perturbations that are not too strong, naturally. In contrast to fiat-band models, the defined 
nearly-fiat-band models can, in principle, be realized in experiments. 

The stability proof was a real breakthrough, despite the fact that some numerical evidence 
and non-rigorous arguments for the stability of ferromagnetism in such models had been 
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earlier obtained [HI [71 [8]. The only comparable result, but referring to an unsaturated 
ferromagnetism (sometimes referred to as ferrimagnetism) , was obtained by E. Lieb [9j. 
While it is hard to overestimate the importance of the above mentioned results for our 
understanding of ferromagnetism, one should bear in mind that the described phenomenon 
is rather fragile, in the sense that it is very sensitive to details of the underlying lattice and 
matrix tij. In the space of the parameters of Hamiltonian ([T]), even general (not necessarily 
ferromagnetic) nearly-flat-band systems are of "measure zero". 

There are many theoretical examples of nearly-fiat-band systems, yet we do not know 
of any measurements performed on real nearly-fiat-band systems. There have been a few 
proposals of experimental realizations of nearly-fiat-band systems, as atomic quantum wires 
[10] , quantum-dot super-lattices [llj or organic polymers [I2]. However, the most promising 
seems to be a realization as an ultracold gas of atoms in an optical lattice [13]. Due to a 
very good control of system parameters in the latter case, such a realization would open new 
possibilities of investigating the mechanism of nearly-fiat-band ferromagnetism, inaccessible 
in other experimental realizations and beyond the scope of present-day theoretical methods. 

Most of the models of nearly-flat-band ferromagnets considered in literature, in particular 
those studied in this paper, represent insulators or semiconductors [2]. Some physicists might 
consider this feature as a drawback of those models. However, recent investigations show 
that ferromagnetic insulators or semiconductors are very promising materials for vigorously 
developing spintronics [14J. Therefore, the mentioned feature constitutes another good reason 
to study those nearly-flat-band electron systems that are non-metallic. 

From the perspective of experimental realizations of nearly-flat-band ferromagnetic sys- 
tems, it is important to know not only the flat-band conditions on tij, which guarantee that 
the lowest band is flat, but also to have an estimate of the threshold value of the Hubbard 
repulsion, Uth, for given widths of the broadened flat band, which are related to the strength 
of a perturbation. It might happen that, for physically realizable parameters tij, the cor- 
responding threshold value Uth is beyond a range of accessible values. While the flat-band 
conditions are known for numerous lattices in D = 1, 2, 3 dimensions, the estimates of Uth 
are scarse. To the best of our knowledge, there are just two papers in literature where this 
question was addressed. The flrst one is a rough estimate by Tasaki [15j, based on his theo- 
rem, with no attempts to optimize it. The second one is by Ichimura et al [16j, and amounts 
to calculating the threshold value of U for a few small systems and then extrapolating the 
results to larger systems. The problem is (see the sections that follow) that a small-system 
threshold value of U, from now on denoted f//^*', depends signiflcantly on the perturbation 
used to broaden the underlying flat band and on boundary conditions applied to the system. 
We show that not only the values of U^f^^', calculated in (TB], but also their least square 
extrapolation to inflnite systems are much above the value of Uth- 

In this paper we extend a little the method of estimating Uth that is based on Tasaki 
theorem. We calculate also t//^*' for small Hubbard systems with open and periodic-boundary 
conditions and flnd, in particular, that contrary to what is suggested in |16j it is not the 
open-boundary condition but the periodic one that is appropriate for the task of estimating 
Uth- 

In order to minimize the computational burden, we have performed calculations for two 
periodic, one dimensional, nearly-flat-band ferromagnetic Hubbard models, whose lattice 
consists of a small number of sublattices. The flrst one, known in literature as A-chain or 
saw-tooth chain (two sublattices), belongs to the Tasaki class of models and a version of this 
system has been analized by Tasaki (see [S] and references there). The second one (three 
sublattices), proposed in [T7] as a model of quantum wires, is off the Tasaki class; we name 
it the sparse A-chain. An attempt at estimating its Uth has been made in [T6j. In both 
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mentioned models, the underlying lattice is a kind of a decorated one-dimensional lattice, 
which is often visualized as a quasi-one-dimensional lattice. 

The plan of the paper is as follows. In Section 2 we give some details of the two methods 
of estimating Uth, used in our paper. Then, in Section 3, we present our results for the 
A-chain, and in Section 4 - the results for the sparse A-chain. A resume of our results can 
be found in Section 5. Finally, in Appendix, we provide a few tables with numerical values 
of some data calculated by us, some of which are depicted in diagrams of Sections 3 and 4. 

2 The methods of estimating Uth 
2.1 The Tasaki method 

The Tasaki method of estimating Uth is based on his proof that a class of nearly-flat-band 
Hubbard systems has the ground state, which is a saturated ferromagnet [5]. This method 
enables one to calculate a value of the Hubbard repulsion, denoted here Uth, above which 
the ground state is a saturated ferromagnet, that is f/^/i is a rigorous upper bound for Uth oi 
sufficiently large lattices. In this subsection, we consider flat-band Hubbard models with the 
periodic-boundary condition, whose lowest band is fiat and is separated from higher bands 
by a gap. The Hamiltonians of such models we denote by Hfh{0), where zero stands for zero 
width of the fiat band. It is convenient to choose the zero of the energy scale at the energy 
of the fiat band; as a result, the Hamiltonian Hff,{0) becomes positive semi-definite. 

In many models, including the two studied in this paper, all the constructions, presented 
below, are feasible. Let the underlying lattice A consist of n sublattices, denoted A, B, C, 
. . . , and let each sublattice consist of JV sites (i.e. JV is the number of elementary cells). Let 
T be the minimal set of lattice translations r such that the set of translated sites, {rilreT; 
coincides with the sublattice to which i belongs. T consists of the multiples, r^, of the 
elementary translation, Tq, m = 0, . . . , A/" — 1, where Tq stands for the identity. 

One can construct a basis of the fiat-band eigensubspace that consists of A^ localized, over- 
lapping (so in general nonorthogonal) states, with supports {Ta}reT, where a is a bounded 
and connected subset of the lattice; the union of these supports coincides with the whole 
lattice. We label the states of this basis by the sites of sublattice A and denote the cre- 
ation operator of an electron in a state belonging to the localized basis of the fiat-band 
eigensubspace, with spin projection a, hy a] ^, i & A. 

After that, one can construct a basis of the orthogonal complement of the fiat-band 
eigensubspace that consists of {n — l)J\f localized states. A/" states with supports {r/3}T-er5 
labeled by the sites of sublattice B, M states with supports {T7}T-er; labeled by the sites of 
sublattice C, . . . , where /3, 7, . . . , are some bounded and connected subsets of the lattice. 
The creation operators of electrons in those states with spin projection a are denoted by 
b]^, i E B, c] ^, i G C,. . . , respectively. By means of the defined creation operators and their 
adjoints, the Hamiltonian Hfb{0) can be written in a manifestly positive semi-definite form: 

after adjusting the amplitudes of the localized states defining the creation and annihilation 
operators. 

Apparently, the two states Hie^ '^lo-IO)) differing by a, belong to the ground-state sub- 
space of the system ^ for N^, = J\f. Their total spin is maximal, Stot = A/'/2, and the value 
of the total-spin projection is for one of them maximal, Af/2, and for the other - minimal. 
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—Mjl. By applying to the state with the maximal total-spin projection the operator S^^^^ up 
to A/" — 1 times, we obtain IStot — 1 states with the total spin Stot = A/'/2 and the total-spin 
projection between Stot and —Stot- All the mentioned states are linearly independent, and 
they constitute a basis of a {2Stot + l)-dimensional subspace of A/'-electron state space, that 
we name the ferromagnetic subspace and denote M ferro- The proof that the ground state 
is a saturated ferromagnet amounts to proving that A4 ferro is the ground-state subspace. 
Mielke and Tasaki succeeded in constructing specific examples of Hamiltonians Hfb{0) (that 
is specific lattices and suitable matrices tij for Hamiltonian ([T| (see [2] and references quoted 
therein), such that the subspace A4 ferro is the ground-state subspace for any U > 0. In their 
proof, a connectivity of the underlying lattice and some overlap properties of the supports 
of the flat-band eigenstates play a key role. 

There are many ways of perturbing fiat-band systems, whose ground state is a saturated 
ferromagnet for any f/ > 0, to get a nearly-flat-band system, whose ground state is a 
saturated ferromagnet only for sufficiently large U. However, in Tasaki method a specific 
perturbation is used [5j. The Hamiltonian of a nearly-fiat-band system,ifjfc(s), obtained by 
perturbing in Tasaki way the considered above fiat-band system with Hamiltonian Hff,{0), 
reads: 

Hfb{s) = -sY^ 4,a<^c,a + XI H^^i.'^ + 5Z ^-^Ci,^ + ■ ■ ■ + ^ ^i^^rii,.^, (3) 

ie^.cr ieBjCr iGC,a iGA,a 

where s > is a parameter measuring the strength of the Tasaki perturbation that broadens 
the fiat band; in the models studied in the sequel s will be related to the width of the 
broadened fiat band. The advantage of the Tasaki perturbation, with respect to other ones, 
is that the ferromagnetic ground-state subspace of the fiat-band case (s = 0), A4 ferro, is an 
exact eigensubspace of Hfb{s), for any s. The corresponding eigenvalue - the ferromagnet's 
energy - amounts to Eferro = —se^/, where e is given by the anticommutation relation, 
{a\^,ai^a} = £, for i & A. Remarkably, the energy Eferro is the lowest eigenvalue of Hfi,{s) 
and Ai ferro is the ground-state subspace, provided that s is small enough and U is large 
enough. This statement has been proved by Tasaki jSj [15] for models that belong to the, 
defined by him, class (nowadays named the Tasaki class) of Hubbard models. Below, we 
outline the proof and simultaneously present the Tasaki method of estimating Uth- 

Choose a bounded and connected subset, S, of the underlying lattice that contains ma 
supports of the fiat-band localized eigenstates: a, Toa,. . ., Tq"^" ^^a; these supports are 
labeled by sites in a subset As^a of sublattice A. Moreover, we require that S contains all 
the supports r/3, T'y, . . . , that have non-void intersection with the union of the supports 
{r™a}m=o,...,ma-i- These supports are labeled by the sites in a subset Bg^p of B, a subset 
Cs^'y of C, . . . , respectively. 

Define an operator Ps - the so called potential, whose support is S, such that Hfb{s) = 
'YliT&T -^tS- following form: 



ma , — ' ' ^ — ^ ' 2 

+ X C{ma)c\^Ci^a + \uCu{,rna)ni^arii-a 



(4) 



+ ... 



where the model-dependent coefficients B{ma), C{ma), ■ ■ ■ , are unique. On the other hand 
there are numerous choices of coefficients Bu{ma), Cu{ma), ■■■■ One can fix the latter 
coefficients by requiring that there is a Hubbard repulsion term at each site of S. We 
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adopt this convention in our calculations carried out for specific models. The potential Ps, 
expressed in terms of operators Cj^, Ci^a, reads: 

Ps = J2 — C'lC',. + I J2 -n^,.n,,-., (5) 
^-^ rrii j ' 2 ^-^ rrii 

i,jS:S,cr '-^ i,a 

where m^j- and mj, i,j G iS, are suitable natural numbers. We note that Adjerro is an 
eigensubspace of Ps to the eigenvalue Esjerro = —se. 

Suppose that for some value of the perturbation, sq, and the Hubbard repulsion, Uq, the 
following conditions are fulfilled: 

(i) Esjerro IS the lowcst eigenvalue of Ps; let Ais be the corresponding eigensubspace, 

(ii) any state |$) G M.Si with Nf, electrons, is of the form 

I*) = E n (6) 

for some states |$(ti,...,o-™.^) with A^^e — rric, electrons; this condition implies that 

i<f) = 4J<^>..) + «Ui<^'.-.), (7) 

for any site i G As,a, where |$i,±o-) are some states with A^^e — 1 electrons, 

(iii) rii^o-ni-al^) = 0, for every site i G 5 and any |$) G Ms- 
Then, Aiferro IS the grouud-statc subspace of Hfb{so) for the specified value Uq of the 

Hubbard repulsion. It is worth to note here that this statement holds true as well for any 
U > Uq, since for such values of U, the [/-independent subspace Ais continues to be the 
ground-state subspace of Ps- Therefore, Uq is a rigorous upper bound for the threshold value 
of U for any system whose lattice A contains S, in particular for a macroscopic (or infinite) 
system. 

If a sequence of upper bounds Uth, obtained for a sequence of potentials whose support 
increases indefinitely, is convergent, then its limit coincides with Uth of a macroscopic system. 
The conditions (i), (ii) and (iii) have been proved to hold true, using analytic arguments, 
in the case of Tasaki class of models [H], for a specific potential; this was achieved in an 
asymptotic region of parameters (where, in particular, U/s is sufficiently large). 

We note that it is enough to investigate the potential Ps in the Fock space corresponding 
to the system whose support is S, and then Ps can be represented by a finite matrix. If 
the support S is small enough, then a numerical exact diagonalization procedure is feasible, 
and the above conditions can be verified numerically. By carrying out such a procedure 
for a given value of perturbation, Sq, and a decreasing sequence of values of Uq, an optimal 
estimate of Uth can be determined. Eventually, one ends up with a computer assisted proof 
of ferromagnetism for given sq and any U > Uth- 

We follow the route, outlined above, i.e. calculate possibly good estimate of Uth for a 
sequence of potentials Ps with increasing supports, first in Section 3, where a model from the 
Tasaki class (A-chain) is studied, and then in Section 4, where a model off the Tasaki class 
(sparse A-chain) is considered. In particular, in the latter case the carried out calculations 
amount to a computer assisted proof that the ground state of a nearly-flat-band sparse 
A-chain is a saturated ferromagnet. 

Actually, the matrices that appear in our calculations are huge and sparse; their size 
reaches the order 10^. To diagonalize such matrices we used the C library Primme of Andreas 
Stathopoulos and James R. McCombs [H], |19].|20]. 
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2.2 Small-system estimates of Uth 



A straightforward way of estimating Uth is to calculate it for systems on small lattices, whose 
number of sites, A^^ constitute an increasing sequence, terminating usually at the highest 
level for which a numerical exact-diagonalization procedure of the underlying Hamiltonian 
is feasible for our computer facilities. We specify a small piece of the underlying lattice, the 
parameters of Hamiltonian ([T|, a boundary condition and the electron number A^^e, and carry 
out an exact numerical diagonalization. The latter procedure provides us, in particular, with 
the lowest eigenvalue and a basis of its eigensubspace. In order to qualify the ground-state 
subspace of our small system as a saturated ferromagnet, we verify if 

(a) the number of elements of this basis is + 1, 

(b) all the elements of this basis are eigenvectors of the square of the total-spin operator to 
the same eigenvalue Ne/2[{Ne/2) + 1] (in units where h = 1), that is the total spin of the 
ground-state subspace is Stot = 

Repeating the above verification procedure for a decreasing sequence of values of U, we 
determine a value of U, denoted U^f^^', that separates the values of U for which the system 
is a saturated ferromagnet from those for which it is not. For a given flat-band model, U^ff' 
depends on A^s, a kind of perturbation and boundary conditions used. Finally, having a 
sequence of the values of U^f^^', obtained for systems of different size, we try to guess how 
U^f^^' scales with Ng, in order to predict the corresponding value, Uth, for a macroscopic 
system. 

We end up this subsection with a comment on the choice of the electron number Ng in 
small-system calculations. Most of the proofs of fiat-band ferromagnetism and the Tasaki 
proof of nearly-fiat-band ferromagnetism can be carried out only for particular A^^e, such that 
exactly half of a band in the single-particle spectrum can be filled. For periodic boundary 
conditions such a amounts to Ng/n. However, for open boundary condition this definition 
looses its meaning, since the bands are not well defined and, moreover, Ng might not be 
divisible by n. In the latter case, we choose one of the two integers that are the closest 
to Ng/n. As will be demonstrated in Sections 3 and 4, the results of calculations depend 
significantly on A^^, whether it is divisible by n, and on the chosen integer approximation to 
half-integer Ng/n. 




Figure 1: The nearly-fiat-band A-chain. The expressions across the bonds of quasi ID 
lattice are hopping intensities, those by the sites are external potentials adjusted to fulfill 
the fiat-band condition when s = 0. 
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3 The case of A-chain 



The first system for wliicli we consider tlie problem of estimating the threshold value Uth of 
a macroscopic system is the A-chain. We shall argue that an extrapolation of small-system 
estimates obtained for periodic-boundary condition is a reasonable approximation to Uth of 
a macroscopic system. The A-chain is a periodic Hubbard system whose quasi ID lattice 
and the parameters of the one-body part of Hamiltonian ([T| are depicted in Fig. [T| for the 
system perturbed in Tasaki way. 




Figure 2: The localized fiat-band eigensubspace basis - a state a|^|0), i G A, and the 
localized basis of the orthogonal complement of the flat-band eigensubspace - a state bl^\0), 
i ^ B. Open circles - the support a - "V- valley", solid circles - the support /3. The 
expressions by the sites are the components of those states. 





Figure 3: A potential whose support consists of 5 sites -one-V-valley potential. The ex- 
pressions across the bonds are the hopping intensities, those by the sites and independent of 
U are external potentials, and those proportional to U are the coefficients in the Hubbard 
repulsive interaction . 



The underlying lattice consists of two sublattices: one formed by the sites in the bases 
of the triangles made by continuous bonds, denoted A, and the other - by the tops of 
these triangles, denoted B. A bond between two sites is represented by a continuous line. 
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if the corresponding hopping intensity is nonzero for any s > 0, and by a dashed hne, if it 
vanishes in the unperturbed system. The hopping intensities between neighboring sites of an 
elementary cell and between sites of two such cells are independent parameters. We choose 
one of them as the energy unit. In the unperturbed case (s = 0), the hopping intensities are 
t+, t_,l, and the on-site external potentials are adjusted to fulfill the flat-band condition: 
for i G A, ti^i = ^ + for i & B, ti^i = t^t-, and to set the fiat-band energy to zero. The 
fiat-band is the lowest one, if > 0. The width of the broadened flat band is not greater 

than 4s and the gap above it is not smaller than (1 + s) (j^ + + ~ 2 j > 0. 

The fiat-band eigensubspace of the defined system is generated by linearly independent 
(nonorthogonal in general) localized states a|o.|0) - "V-valley", for i E A. The support of 
such a state, a, and the corresponding components are shown in Fig. [2j A localized basis of 
the orthogonal complement of the fiat-band eigensubspace consists of states denoted 6jg.|0), 
i & B. The support of such a state, /3, and the corresponding components are also shown in 
Fig.|2} 




y 1 V 1 V 1 V 

3 3 3 

*+ *+ I _ „ *+ I *" _ „ *+*- *- 

3t_ 3t_ 3t+ ''2 3t_ 3t+ ^2 3t+ 

K U. V. H 

4 4 4 4 



Figure 4: A potential whose support consists of 7 sites - two- V-valley potential. The 
expressions across the bonds are the hopping intensities, those by the sites and independent 
of U are external potentials, and those proportional to U give the coefficients in the Hubbard 
repulsive interaction. 



To apply the Tasaki method, we define the potentials Ps'- 

U 



S I 



(8) 



U 



2(m„ + 2) 



t.CT ' ''2,— (T 



Specifically, the potentials (|8]) containing one- V-valley and two- V- valleys, are given, in the 
form ([5]), in Figs. 3]|4 
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Figure 5: The case of a nearly-flat-band A-chain, with t+ = 1.0 and t_ = 1.5, perturbed in 
Tasaki way, with s = 0.01. Uth versus the inverse number of sites, 1/Ns, in the supports S 
of the potentials Ps, defined in dsj), with A^^ = 5, 7, 9, 11. 



Choosing a specific A-chain and calculating Uth for the defined potentials, whose supports 
consist of A^^ = 5, 7, 9, 11 sites, gives us a decreasing with A'^ sequence, which might saturate 
at a nonzero value. Fig. [H] The smallest Uth is our best estimate of Uth obtained by Tasaki 
method. 

In Fig. [6] we give the values of U^f^^' obtained for small systems with open and periodic- 
boundary conditions; for comparison we include also in this figure the values of Uth obtained 
by Tasaki method. Clearly, the values of f//^*' obtained for open-boundary condition are 
much above the rigorous upper bounds provided by Tasaki method. The values of U^f^'^' 
obtained for periodic-boundary condition seem to oscillate around a constant and are a little 
below the values of Uth- 

The behaviour of the latter two kinds of estimates of Uth is more clearly seen in Fig. [7| It is 
tempting to suggest that the best (say, in the least-square fit sense) constant approximating 
the values of Uth obtained for periodic-boundary condition is a good approximation to Uth, 
lower than the lowest value of Uth, for the considered sizes of small systems and supports of 
the potentials. 
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TasaJd method, Ns=5,7,9,ll. 

Small systems with periodic b.c, 
Ns=6,8,10,12, Ne=3,4,5,6, respectlYely. 

Small systems with open b.c, 
Ns=5,7,9,ll, Ne=3, 4,5,6, respectively. 

Small systems with open b.c, 
Ns=6,8,10,12, Ne=3,4,5,6, respectively. 



0.14 

1/N 



Figure 6: The case of a nearly-flat-band A-chain, with t+ = 1.0, t_ = 1.5, perturbed in 
Tasaki way, with s = 0.01. The threshold value of U versus the inverse number of sites 
1/A^, i.e. f//^"*' versus l/Ng - for small-system data, and Uth versus l/Ng - obtained by 
Tasaki method with the potentials defined in (|8|. Numerical values of the data for periodic 
boundary condition (periodic b.c.) are given in TabjlJ 



0.14 

1/N 



Figure 7: A part of data presented in Fig. ([6]), displayed with a higher resolution. The 
case of a nearly-fiat-band A-chain, with = 1.0, t_ = 1.5, perturbed in Tasaki way, with 
s = 0.01. The threshold value of U versus the inverse number of sites 1/A^: solid squares - 
U^f^'^' versus l/Ng for small-systems with periodic boundary condition, bullets - Uth versus 
1/Ns obtained by Tasaki method with the potentials defined in Q. Numerical values of the 
data for periodic boundary condition (periodic b.c.) are given in Tabjl] 
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4 The case of sparse A-chain 



The second system, we consider, is a modification of tlie A-cliain, wliicli we name the sparse 
A-chain. This is a periodic Hubbard system whose quasi ID lattice and the parameters of 
the one-body part of Hamiltonian ([T]), with Tasaki perturbation, are depicted in Fig. [s] 

An attempt to estimate Uth of a macroscopic nearly-flat-band sparse A-chain was made 
in |16) . The underlying lattice consists of three sublattices: the sublattice A that consists of 
the left sites of the bases of the triangles made by continuous bonds, the sublattice B that 
consists of the the tops of those triangles, and the sublattice C that consists of the right sites 
of the bases of those triangles. A bond between two sites is represented by a continuous line, 
if the corresponding hopping intensity is nonzero for any s > 0, and by a dashed line, if it 
vanishes in the unperturbed system. There are three independent hopping intensities. The 
hopping intensities between neighboring sites of the sublattice C and the sublattices A and 
B are the same, t, in an unperturbed system. Then, the hopping intensity between two sites 
belonging to the base of a triangle made by continuous bonds is chosen as the energy unit. 




m + st 



m + sf^ 



m + 1 — si^ 



TO + 1 — s< 



TO + 1 — st 



TO, + 1 — st^ 



Figure 8: The nearly-fiat-band sparse A-chain, for Tasaki perturbation. The expressions 
across the bonds of quasi ID lattice are hopping intensities, those by the sites are external 
potentials adjusted to fulfill the fiat-band condition, when s = 0. 



Finally, the hopping intensity between two neighboring sites in the bases of neighboring 
triangles is m, in an unperturbed system. The hoping intensities, t, m, are independent 
parameters. On the other hand, the on-site external potentials are adjusted to fulfill the 
fiat-band condition and to set the flat-band energy to zero: the potential at the sites of the 
sublattice C is and that at the sites of the sublattices A and B \s m+ 1 (see Fig. [s]). We 
shall see below that the fiat-band is the lowest one, if m > (t can be of any sign). In order 
to make our model the counterpart of that in |16| . we have to set an extra relation between 
t, m: t = y/m for m > 0. The calculations have been carried out for m = 0.8 (the value 
chosen by Ichimura et al [IS]), and s = 0.070. For this particular sparse- A-chain, the width 
of the broadened fiat band (Tasaki perturbation) is 0.136, which is 2s approximately, and 
the gap above the broadened fiat band is 0.342. 

The support a - of a localized element of the flat-band eigensubspace basis, al^\0), i ^ A 
and its components are shown in Fig. [9j Such a support can be depicted as a valley between 
two neighboring triangles made by continuous bonds - a "W- valley". The supports, /3 and 
7, and components of elements of a localized basis of the orthogonal complement of the 
flat-band eigensubspace, bl^\0), i G B, and c|o.|0), i G C, are also shown in Fig. 9 
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t 




t -t 



Figure 9: Localized flat-band basis - the states aj^|0), i G A, and localized basis of the 

orthogonal complement of flat-band eigensubspace - the states 6j^|0), i & B, and c|g.|0), 
i & C. Open circles - the support of a localized flat-band eigenstate ("W-valley"), open 
squares - the support of 6|^|0), i G i3, open rhombs - the support of c]^^|0), i G C. The 
expressions by the sites stand for the components of those states. 



By means of the defined above fermion creation operators and the corresponding anni- 
hilation operators, the Hamiltonian of the flat-band sparse A-chain, can be written in the 
positive semi-definite form ([2]). This proves that m > is a sufficient condition for the flat 
band to be the lowest one. 

The potentials Ps that we use when applying Tasaki method are defined as follows: 

Ps = -— + — + 

irta , rria ^-^ ma + 1 

i&As,a,cr i&Bs.i3,cr i&Cs,^,^ 



u ^ 

H > ni fjUi 



In the form ([s]) they are shown in Figs. lOpT , for A^^ = 6,9. One can check easily that if 



the conditions (i), (ii) and (iii) of the Tasaki method are satisfied for these potentials, then 
the Tasaki proof of saturated ferromagnetism in the ground-state [5] can be carried out. 

Up to now, we have dealt with only one way of perturbing a flat-band system to obtain a 
nearly- flat-band one, i.e. the Tasaki perturbation defined in ([s]). An alternative way consists 
in modifying one or more arbitrary chosen hoping intensities and/or external potentials, 
while keeping the remaining ones intact, so that the flat-band condition is no longer valid. 
This is the way adopted by Ichimura et al in [16j to perturb the above defined flat-band 
sparse A-chain. 

Apparently, different perturbations of the same flat-band system result in different nearly- 
flat-band systems. In order to compare the threshold values of U calculated for such a 
different systems we have to make the systems itself comparable. We have chosen to consider 
different kinds of perturbations as equivalent, or having the same strength, if the widths of 
the broadened flat bands in two nearly-flat-band systems, obtained from a given flat-band 
one, are the same. To motivate our choice, we note that according to good qualitative 
arguments [8j the way of perturbing a flat-band system is not relevant for the existence 
of the phenomenon of nearly-flat-band ferromagnetism; whatever the perturbation, there is 
Uth > 0, which increases with the strength of perturbation. 
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Figure 10: A potential whose support consists of 6 sites -one-W-valley potential. The 
expressions across the bonds are the hopping intensities, those by the sites, and independent 
of U, are external potentials, and those proportional to U are the coefficients in the Hubbard 
repulsive interaction. 



Moreover, we expect that for sufficiently large electron systems, with sufficiently narrow 
nearly-fiat-bands, the threshold value of U is an increasing function of the width of the 
nearly-flat-band only. In other words, Uth is an intrinsic property of nearly-flat-band elec- 
tron systems, independent of boundary conditions or a kind of weak perturbation applied 
to the underlying fiat-band system; the only relevant physical parameter, characterizing a 
sufficiently narrow nearly-fiat-band, is its width. 




Figure 11: A potential whose support consists of 9 sites - two- W- valley potential. The 
expressions across the bonds are the hopping intensities, those by the sites, and independent 
of U, are external potentials, and those proportional to U are the coefficients in the Hubbard 
repulsive interaction. 



In |16| . the fiat-band sparse A-chain is perturbed by replacing the hoping intensity t = 
\fm by t = \fm -|- 5, for some 5 7^ (the Ichimura et al perturbation). Now, we are 
ready to compare estimates of IJth for two nearly-fiat-band sparse A-chains, one obtained by 
perturbing the fiat-band sparse A-chain in Tasaki way and the other - in Ichimura et al way; 
by adjusting the parameter s of Tasaki perturbation we can make equal the widths of the 
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broadened flat bands in the both cases. Specifically, for m = 0.8 and t = 1 the width of the 
broadened flat band (Ichimura et al perturbation) is 0.136, and with Tasaki perturbation 
the same width is attained for s = 0.070. The gap above the broadened flat band is 0.210. 

In Fig. 12, we depict a collection of the results of our calculations of threshold values 
of U, obtained for perturbed flat-band sparse A-chains. The rigorous upper bounds, Uth 
versus l/N^ (see Tab. |4]for numerical values), obtained by Tasaki method, constitute a set 
of reference data. This is a strictly increasing sequence, with a tendency to saturate for small 
values of l/Ns] very much like in the case of A-chain. 

The threshold values U^f^^' for periodic boundary condition versus l/Ng seem to oscillate 
around a constant value, irrespectively of whether the system is perturbed in Ichimura et al 
way or Tasaki way. The numerical values of these data can be found in Tabs. [2| |3j In both 
cases, of Ichimura et al perturbation and Tasaki one, the data are below the rigorous upper 
bounds. 



u3 



0.08 



0.10 



0.12 

1/N 



0.14 



0.16 



Figure 12: The case of a nearly-flat-band sparse A-chain. The hoping intensities of the 
unperturbed flat-band system are: m = 0.80, t = ^/(L80. For the system perturbed in 
Ichimura et al way, with m = 0.80 and t = 1.0, U^f^'^' versus l/Ng-. open triangles - open- 
boundary condition with A^^ = 6,9,12,15 and N^, = 2,3,4,5, respectively; solid triangles 
- periodic-boundary condition with A^^ = 9,12,15 and Ng = 3,4,5, respectively. For the 
system perturbed in Tasaki way, with periodic boundary condition, m = 0.80, t = v^0.80 and 
s = 0.07 (for such a value of s the widths of the flat-band broadened by Tasaki and Ichimura 



et al perturbations are approximately the same) 



TTS.S. 



versus l/Ng, with A^^ = 9, 12, 15 and 



iVg = 3, 4, 5, respectively, - solid squares. Uth versus 1/Ns, with s = 0.07 and the underlying 
potentials defined in Figs. TOpT obtained for A'^ = 6, 9, 12, 15 - bullets. For the numerical 
values of depicted data see Tables. |2} [3| |4j 



In contrast, the small-system threshold values of U obtained for open-boundary condition 



are definitely to high. In the case of Ichimura et al perturbation, the depicted in Fig. [12 
data (which reproduce the results of jl6]) seem to scale linearly with 1/A^s. However, except 
the case of 6 sites, their values are greater than the corresponding rigorous upper bounds 
provided by Tasaki method. Moreover, also the value of a least-square linear extrapolation to 
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1/Ns — )■ limit, is above the rigorous upper bound obtained from a potential whose support 
consists of 15 sites. The Ichimura et al choice of the hopping parameters results in a nearly- 
flat-band width of the order 10~^ of the energy unit. We have repeated the calculations 
for another values of hopping parameters that give the nearly-fiat-band width of the order 
10~^. As expected, the calculated values of are smaller, they also seem to scale linearly 
with 1/Ns, but the slope of the least-square linear fit is smaller than in the previous case, 
so that the value of a least-square linear extrapolation to l/Ng — limit is higher than in 
the previous case. In the case of Tasaki perturbation (open boundary condition) the values 
of U^f^^' are even higher than in the case of Ichimura et al perturbation (see Tab. [s]). 

Apparently, the linear scaling, observed in the case of the fiat-band sparse A-chain with 
Ichimura et al perturbation, is accidental, and the predicted value of Uth is not reliable. 
Taking into account our results for the A-chain, one can claim that small-system and open- 
boundary condition threshold values of U are irrelevant for our task of estimating Uth of 
nearly-fiat-band systems. Definitely, the threshold values of U obtained for periodic bound- 
ary conditions constitute a better basis for deriving reliable estimates of Uth- 

5 Summary 

Nearly-flat-band ferromagnets constitute a speciflc class of strongly-correlated electrons sys- 
tems. Theoretically, they can be thought of as weak perturbations of the underlying un- 
physical flat-band systems. One of their characteristic features, which is a result of the 
competition between the kinetic energy and the Coulomb interaction energy, is the thresh- 
old value Uth of the screened Coulomb repulsion, above which a nearly-flat-band electron 
system becomes ferromagnetic. This quantity is well deflned for sufficiently large systems, 
presumably not necessarily as large as macroscopic ones. It is desirable to have a reasonable 
estimate of that intrinsic property of nearly-flat-band ferromagnets. 

In this paper, we have made attempts at estimating Uth in two nearly-flat-band systems, 
modeled by Hubbard Hamiltonians ([T|, the nearly-flat-band A-chain and the nearly-flat- 
band sparse A-chain, deflned in Section 3 and Section 4, respectively. This choice has been 
convenient, since it enabled us to avoid large volume computer calculations. Concerning the 
electric conductivity both models can be classifled as insulators or semiconductors. 

We have considered two methods of estimating Uth, described in detail in Section 2. One 
of the methods, the Tasaki method, provides us with rigorous upper bounds for true Uth, 
independent of the size of the system and boundary conditions. We have shown how to 
obtain the best bounds, whose quality is limited only by available computer facilities, and 
obtained those bounds for both models. In the case of the sparse A-chain, this resulted 
in a computer assisted proof of nearly-flat-band ferromagnetism (the existence of flat-band 
ferromagnetism was demonstrated in |16j). 

The other method, studied by us, is the method of small-system estimates of Uth', this 
method amounts to what a physicist would typically do when faced with such a problem. 
It provides us with boundary condition and perturbation dependent threshold values U^^^\ 
Then, the problem we have to deal with is how to extract reliable estimates of Uth, having 
a set of values of U^f^'^' obtained for systems of different sizes, open- or periodic-boundary 
conditions, and various perturbations, like Tasaki or Ichimura et al perturbations. Our main 
conclusion is that when applying the method of small-system estimates one should resort to 
periodic-boundary. Concerning the perturbation, the Tasaki perturbation is preferable. 

Finally, concerning the question of the range of values taken by Uth, let us note that, 
in the two models considered, we have chosen one of the hopping intensities as the energy 
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unit. The remaining independent hopping intensities have been chosen to differ from this 
unit by not more than a few ten per cent. Then, the widths of the broadened fiat-bands do 
not exceed 14 percent of the energy unit. The resulting estimates of Uth are either smaller or 
almost equal to the energy unit. However, for the widths of the broadened flat bands that 
are about a few per cent of the energy unit, Uth is significantly smaller than the energy unit. 
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8 


10 


12 




3 


4 


5 


6 


Uth 


0.104 


0.109 


0.105 


0.110 



Table 1: U^(^' for A-chains with periodic-boundary condition and Tasalci perturbation; 
t+ = 1.0, t_ = 1.5, s = 0.01. These data are depicted in Figsjoj [T] (solid squares). The 
least-square estimate of a constant is 0.107. 





9 


12 


15 




3 


4 


5 


Uth 


0.592 


0.602 


0.500 



Table 2: f//^*' for nearly- flat-band sparse A-chains with periodic-boundary condition and 
Ichimura et al perturbation; m = 0.8, t = 1.0. These data are depicted in Fig 12 (solid 
triangles). The least-square estimate of a constant is 0.565. 



Ns 


9 


12 


15 


Ne 


3 


4 


5 


Uth 


0.994 


1.110 


1.094 



Table 3: f//^*' for nearly- fiat-band sparse A-chains with periodic-boundary condition and 
Tasaki perturbation; m = 0.8, t = a/0.8, s = 0.07. These data are depicted in Fig, 12 (solid 
squares). The least-square estimate of a constant is 1.066. 
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Ns 


6 


9 


12 


15 




2 


3 


4 


5 


Uth 


5.900 


3.147 


2.498 


2.225 



Table 4: Uth obtained from nearly-flat-band sparse A-chains; m = 0.8, t = a/0.8, s = 0.07. 



These data are depicted in Fig, 12 (bullets 
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6 


9 


9 


12 


12 


15 


15 




2 


3 


3 


4 


4 


5 


5 


6 


Uth 


22.030 


4.613 


8.986 


7.71 


8.838 


10.195 


9.290 


11.016 



Table 5: U^f^^' for small nearly-flat-band sparse A-chains with open-boundary condition and 
Tasaki perturbation; m = 0.8, t = \/0.8, s = 0.07. 
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